Model Overview

We estimate the age of a site by calculating the years since the last fire. We then fit a curve to model the recovery of vegetation (measured using NDVI) as a function of it’s age. An additional level models the parameters of the negative exponential curve as a function of environmental variables. This means that sites with similar environmental conditions should have similar recovery curves.

Input data

The model was last fit on 2022-03-17 13:19:55.

This version of the model was fit with 1076106 pixels including data from 2002-04-23 to 2019-12-19.

Workflow

This repository was developed using the Targets framework.

── Attaching packages ─────────────────────────────────────── tidyverse 1.3.1 ── ✔ ggplot2 3.3.5 ✔ purrr 0.3.4 ✔ tibble 3.1.6 ✔ dplyr 1.0.8 ✔ tidyr 1.2.0 ✔ stringr 1.4.0 ✔ readr 2.1.2 ✔ forcats 0.5.1 ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ── ✖ dplyr::filter() masks stats::filter() ✖ dplyr::lag() masks stats::lag()

Attaching package: ‘arrow’

The following object is masked from ‘package:utils’:

timestamp

CmdStan path set to: /home/rstudio/.cmdstanr/cmdstan-2.29.1

Results

Environmental Associations with Ecosystem Recovery

These parameters represent the relationship of the following environmental variables to the recovery trajectory.

Recovery Trajectories

The plot below illustrates some example recovery trajectories. It currently just shows the top 20 cells with the most observations.

Model Performance

Compare estimated vs observed values for all pixels. This is not true validation - these pixels were included in the model fitting.

model_prediction %>% 
#  filter(cellID%in%cells_with_long_records$cellID) %>% 
  ggplot(aes(x=median,y=y_obs)) +
  geom_hex(bins=50)+
  geom_smooth(method = "lm",col="red")+
  geom_abline(color="blue")+
#  facet_wrap(~cellID) +
  scale_fill_viridis_c()+
  labs(x="Estimated NDVI",y="Observed NDVI",
       caption = "Blue line is 1:1, Red line is least squares regression. Count is the number of pixels in that location",
       title = "Estimated vs. Observed NDVI"
      ) +
  theme_bw()
## Warning: Computation failed in `stat_binhex()`:
## `geom_smooth()` using formula 'y ~ x'


## Spatial Predictions

Maps of spatial parameters in the model.  


```{=html}
<div id="htmlwidget-ba5659cc4d7ad7328234" style="width:672px;height:480px;" class="leaflet html-widget"></div>
<script type="application/json" data-for="htmlwidget-ba5659cc4d7ad7328234">{"x":{"options":{"crs":{"crsClass":"L.CRS.EPSG3857","code":null,"proj4def":null,"projectedBounds":null,"options":{}}},"setView":[[-33.998707,18.577],10,[]],"calls":[{"method":"addProviderTiles","args":["Esri.WorldImagery",null,null,{"errorTileUrl":"","noWrap":false,"detectRetina":false}]},{"method":"addRasterImage","args":["data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAC0AAABfCAYAAAB88E5QAAAIOklEQVR4nO2beZAU1RnAf+91z/TMzszeu+OyhGNBd7PgeiDBKyjxiIVGhVUMmFh4g5qKMaUmf5gyZQWrxKQwFY3RlGiquFRWxSMeiEYTDRgRCAFFsoC4suuy7MzOsTN9vJc/KAmJWhVnegIx8/tvprve96vvfdOv+3s9UKbM5yJKMahOLNXkM+iBXmT7bb7HkH4POGf2Un3w50l3vqw/79xCMf0c7OIrH9NkHS78SQVm3iJb2USTl/EzBOCjdOe1KzWuwokGUVKgpaCyP4sy/K9AX8pj9e5HNOz/gQilERo8UyK0BgFzZi/ztUSKltZ77jsgJJQmkHEAsMMmRs4jNJhn6bLZvqa7qPLILThfA7TW5HCtKrQUBBN5bMugtifFkhVzSnJ1KkramnsOL53+NCef4xCdNp/UqybZxgqev+u8ksh+QsHSaufC/VmeINn840vxbnVZ9UhnSWU/oeCalmNuFjo5wOiuZ8WdC2rQ4r/iW6ZMmTJlypQpU6ZMmTJfEPX+At+7R/8JRbcQ1AvXHRLxotCZJ7R66Yb/PfEvBfa9F2vvqStLln3fn/u9lZdrrzeNMbISEQmjU1nktDMQ1Zf6Fsv3/rTRuViYE0dC3oVw6MD3emiF7h9erLV+pegZ8F0agEgYOfUkcD1ERysIA4aTBw6rP/ygKPGSSMsTFgj1+puQySKaWrEjMXQuhWVEcFUeqqJFje/rTsDBiBENqO5eZPdGgi3Hot7eTEV6PeZlS4qu7ZJJA4hoELVxO3L3HsTx7ciqJmBJ8eMWr/b5qK13aOf5d5DVFjJmYVz0cLlLWaZMmS8Le7IPaPfhOSW/ty5qGXfVS1rrNVqt/p4evv3cA7IXLH6mpOIFSw/Zj35KrCE0BtnRwoUnJbj4qsdLJl6QdMreLzRk95FxExCvQ4RNxFPL4bhpvJsQfNwc89f0IApaVtWG23S2/VgUimjvbnZev4Yxt01k+LgTGXaT1Ml6hHVuyZbsL5xpteUODRBauhJDmAzFm6h6bBZi/AQqbMX6vYBp+S56MF9YWv91O8QiePuGyXlpLCMKSD4MZnl5MMVZI5pZ3t1bAtV/UtAUJuefqUMtUYwRUYxzp+NVNpB2BujJpGmPNmMbGsuYfviUB0DVr1eL/rds+lbtQ/9tPYaCkBGl1vJ4Y6CXdX0f++35LxT8EDDi++2I1gkMV9UCWQSSSKCW8VVJknZp15eiptBeNFNrpQmcczwiPgoCIfSmdchTf17Sm/2iVsTgjV0i0BZHNDSz5usrIJtAjD+ShRuWHJ4r4ieoZAa3rplJ71xG5mePMv2pUQC8O/jg4bUiHswzN/fxp457qNq7l/zCuQQtxfH1Ls2RGj/8PpOipadMNdi1LY9a+za7UjY3nb2XWssl4w764feZFC0dv3wc3+kah7rgEgJS80HaJGkbhIwYSq/RuQXna/XcvMPrvTx51q+EqIxg3/JzauZ30RJzOK26kvCiZbC668B5M65/wjdxXy5Ntve8DngalAtODu/Zp/cfcDxkTQS+2sKMRW1YGZtHH5pVdExfenmB3DBDpNDb/kI2ZJKZMR3n3QFy63oZXrMDxk6GSSZCwezvLi86475kWr3+Q01DDSJWDUe0ktc5zMcfo39VP/EbJiKmnE3KSxB9521Ey1GkqyqJBS8qOLYvmVbb96C3/h2n61X0y09iKQO3O8Hvl+Xxtu6GRA/RQB35SaeQqa4GIG2vLDjjvjQgRXUIb1cClXFwNvXApgcw6sJ0zssi6yqg8ghclcf2snjawdUOhjAKjueLtDHjIaHW3qrFUUfjLHkSlbI5OXk1zUaKxPowr0ztx6huIhKoIeXs5ZqbazDzXsHxfGuqi5ZWsGLkrplF+sZOavfl0VIwcusAZ/yiHikMtFZYMsIHzTGckMHM+V0FlYh//elQFN3/PrFoI7l7VnDdtbNpqnAYHdPc9EKWmdcP41omNd1J2gMpnGgQL1hYznzLtE724o2cSDoWJnTrVXwrtZ3J2qExNJb7prtMmDuME5Bk4xGcaBDpKqxkvqBYvt33Og9+W8t4DHFMK+z4EOfUs5DCwEztQ/ftQFTWgpC86IT45aJaDFex6p4LCorvX3mYAnvDHgLJLHLiWKQ0SeZ7qQvXod/bgRraijy5gzM3rqXt9imMjl17aK/TAEa8CmvKKETUIvPbP2NsfoNqqwlvxTK8niQYAnb1IMIWX0nk0LnCW2e+ZVpOv194T1+tdcbGdQSieRwCgZg2GalcRLwV8mn0hrWwaQu67qOCY/m6u+Vu7UcNO1iNJjrVh6yogsYj0cEw3vLFyNZRiBPPRL/3FrLjp4e+PACsW54URm0Ya/IIyNsgA7imgaNyyEvmIo49A4TBa3Vj6ZxX+DLu+46t2RYHT6G37UQ99zhmah9BgggnB33bsBf9htMsCKadgmOU5FFfvXmLJhREd/fgbOnHiEfAlOQ39vNxd4gxP2pDnnJ3wbFL1p9w1YtabluHu2Yj6fVJKsaFwFX7W2lfa4O8jTxhwSG+Tv8b8oUnoH08RjxC1XmViDHNiIZG9M6drAiNZ6iIwizd3viIBuzfvUKgrQ553kzu2pLgjwujeOZo0tEgr93+zYJnuXTSHacTOHoq927uY8ndFiG7niMSScycS6CxoqihS/OSCiDFN4QUBjfUDnD/NbswlCIfDTI4MkaiJlzc2D45fiZCTBMiFOWYzD7qP0oTzNrkLRNPFvf7L6k0gKi/Qoh4K+kpEazBPPUfDlE3kC1qzJJLA4jw+eLojkGcWJBA1iXa5///bktG57yVetYVn95/LFPm/51/AFKYGcGgIviuAAAAAElFTkSuQmCC",[[-33.941641403099,18.3003829808018],[-34.3371886093941,18.5267584324]],1,null,null,"alpha"]},{"method":"addRasterImage","args":["data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAC0AAABfCAYAAAB88E5QAAAHh0lEQVR4nO2bW2wcZxXHf983szve9XrXXntjO04anMglSoIIKSJtmiCKSAr0pgopAkGAoj6ABDwhhKoqV5IgghCIIpUHHqqAkKpWBJJyEaRVUhFBlNLSOJc2F9upL3G83l3be5vb9/EwcQhtitTd2SRF+3uxd2Sf8/eZMzPnMoYmTd4V0Qij6vj3NJ4HloW8a3foPmTYBrOVZ/T1n9XwPv1uP1srZpjGxktPawB7zTry9jiZln5KXj5MF0CIog+PPqOhgqts/jo6gxBx4sZl5tzQT2Y46aHe2K0B/nQpwaGROLOOIGEq8o7ByZxk3a6/hJoidUdandyuAWKmwtMGbRFY3m4zNGfxzJF22g/kOPa7jaFejHUZU0M/1BTLkM3zQv9d9LW6nM5bJKMay1Dcv/grDbk71RfpuRIi0wPdizh01sAutfD46hLrex9riNh5ajauBndolILOLv7gpjgyYbDvni82VOw8NV+IctU2gWkg+74lXhw3GEiFfjtu0qRJkyZNmjRp0qRJkxBQp3bcknal7rmHOrXz/ddnqVM7tDq39/0n/P8CNbxPq3N7Ghb90OcU6s09mmgEXBeUBiGgWEKuCW9O3Zih+uB2jWWB74OUUCwF4ienoacLuXpXXX7Dn8NCIDCbB6XA9YKIm/+ZwKlLP64rdRoiWq7cJljYDW1JMI1AeNUGKwqeD8qvy36om4D/wvdBGsH3tgM6Am2tyI/uqTslGycawLWDFIlG0BNTqHw1FLMNnXKq8z/QGBJmi1CuIO/Zd1Omqk2aNGlyE6juelDP7xUbSV0PF2//l7RIWog7eqE9ifOro8Hx32zR5hf2N+yeXHPtUd378DsiavS0QjaPbIvh/vLzDYt4TaLdpzdrAGdwCn+4ABUbhMS4bzXu8WFmDw4h4pFwlV5HTelh3LsCeXoIla+gPQWpBGRzMDmN0dNKcm0/8hM/uX3SQ72+TQOIDw8gkhb2mQKvfPUE5GfBijLy7Cy0xsJXeh3vPdJSBmVnxMQbmUFG4Py5Com9l7gy5bHhFys58MgrDZB6nYT3/AurtgmsKJQqRD+1EqOzhU2fkzg2uJ5GXZjgge8kGqH1GrUv9M/v1SgNI+OoywXyh6d49ZhEA+1Jg7UnDt0+OX2NqRxcfAs1kUfN2cQzmnhMMj3pkogbIUp8J/W9pHJ6p8bzIZtHF8vosoPsTCI3PtXQYr+uxlau2CoQAuIxxKoBZLoNopEbPnjCpP5uvFxBZ/M4zx6D9jbE2k8CcO7TDzVMeP2NbXcXoiuNOV3CPXqWSCbNzGmHv7/shiDvxtQfaduBSBSZacOfqjD0zaO89IKHW25chtQvumJDNguGxEi34JQlrq3pWGgy/bVNWr2+Teuxn4X6F4RylauT2zXJNvA9/MOvYixdAN3pYHaXSSM6u6AliegM55W3cMZicyUolaDqYHzsTsh0QCQCHSm4kkO/eR4iLejygVAiHoro6sGTMDIOaLicRQ+NBrneGoPOFIVfn4HKLBgRdOm3dQsP5XSNbv6sLs8YLLnfJLJ+OeQKkGqDdCo4C4VZdKGI6FsAjguOg9z485p9hxLpBRviZO6wmRssU95/HFrjQTUoRBDteAzRnoDCHJSDeZ568ds1RzyUAaR59wCplSWwXSb3DxNfsgjGJoKo+n4wVNeAKcH10IUKWLV3NqHVCOrs9zW+D5aFf+Q1zC1fBtOCyiz6ny9DqQpREwyJ3vAA0i4jYo/U5D+8oXrVBttFD55Dl1z0lfPos38Draje/XFYmAnSxfWQs1PgVNDF52pKkdDm03p4HNHfh0inMBf3gFb4/zjDqcdP8KEnehDr1sOdHeQokKuW+EBiMWZ5riZf4Q3VlcL547/AU0RWdCGyecZ/P0uxDM8/NsRnNl8k8dPvYulWLhUVL41nWZqsbY0RXk6f3K5x3GCAHo2gLkzgXiyAEMgWA+34RO/7IPRmyPf1cXQcHl26pSb/4UXaikLEDBZDuRlU0aE67pB8uB+R6UBfmsT+8xmsjZqOi6M8WkejEJ5oQSDYM2BBGtP3aUtEEcsWg5QIz8e4UsI7NYpzZprKzgd1bGttfWRodw858ITAV0G0LQs9U4RoEHXys/hvjDFysIo3Nkd0oAMha8/McLdben6tXA6+Ko2eygdl66olLOvO413IgtK0PHmwZtXh78bnHzIzRfSVaVSugmyNIDLtYBjgKxjox+5dRsx86BZfiPMIEayUM2lEPIZhTaF9/+pxA3q7oFqhpZCt3UWIcq8RRFsFIqUM0kapoA4RAip2Xf8915jd+PInhb741tWKToNx1Y0QwbGp3K2v8m6EWNwLpXIQbccNBANYUcTae8Eu1my7Ma9OAOq187CgM6jqxiYhV0B0d0MqBa6N6Pn6bfBweRuVE5OUD+cwTE16+ybQCmIpdHsfwnPqst2wSLc9dVgkB0zS3/gIjE4AAjoWIV2bIqW6bDdMNEBs6yFBNArJBKIlAXNT6IiFpr7etqGiAeTqnQLXBbMFfBfGB2mbK9dnMyRt/9vJmt0CKfFSGUSsHZ0buxluw0EXn9Nq+EfNtyWbNHkb/wbXBgWo9H39UAAAAABJRU5ErkJggg==",[[-33.941641403099,18.3003829808018],[-34.3371886093941,18.5267584324]],1,null,null,"gamma"]},{"method":"addRasterImage","args":["data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAC0AAABfCAYAAAB88E5QAAAIP0lEQVR4nO2ae4xU1R3HP+e+Zubu7HNm37vsAj5YhCCKkgJpSCuKWBuBEktQWqPBVzFWG9PGtNWmVpumsUkTglGx4m5rUVALViOlFotFNIJUcKUCurLue3d2dndm7sx9nP6xlZBYbZm5t0ozn/9m5s75ffO7v/s753zPhSJFPhURxKAPv9Muu1OC2ojkllnX+B5D8XvAu/d1yFM/X3Hndvlp1+aLr6LX75kU/FqPxpEPTDZtr+fclWN+hgB8LI+fH+iQ3SlIWyoLGx0W1WeoChlIJNWRb/taIr4MtvFwu0zmBFFdMidm01Kaoyel01hiU6pXUNr/EcqUO3wTXnB5eH+5XQK4EkIK9GV0njwa5axynUbLoPRPf/ZVMBQo2vvHzySAqcGX63PseqeEvf0CT8LOboGIXyeUKx7yvXtohfz5Wa2Z5bPqcao3sbBnFQ1LBjmSDHP5lLWBtNKPyXvw7R88IW0PlrfUw1g/KAqifHWgYj8m7/K4svVaEdEkQr1EIF3kWJ+fuooUKVKkSJEiRYoUKVLEJ+55o8N39+i/oWAL4advfj7CC+LHb3TIRzrbzzzh/xds6myXv3k3uOz77lM82tkuBaAIcCSoAvYPCX69aI1vsQIxVx7pbJeGAhkHDHVSdDIR4qOuKLouefHGywuK67upDpOZGLEEVWFIZgWaAlrYOfn76qe3FVQ6gYi+vu0a0RCV2B4IAbURuKDOoaIih5TguoWFLciA/CwmcoCAiCYZtwUzKl0aFoyxYlrhZzCBiYZJv7pUhzrTAyB2SokUQqAu58bD7TKqT3YQTYFV0/0/6SpSpEiRzwk5/Lj09t4V+Nq6oMnl6i3PyG/Om+CqKg052n3y+4yzXUa0KwPryXkvAlZf8+QnMyo9iJqECXFk9OHAMp6X6FU3PD0p6AWbHcd1vLJq6OpBnjhO/7R67j0wyO6eiK9CTyUv0eYKj0xFiFxEZ2qli/jwAGLGTEilqVPquPk8i3Uzg5uyT1v02h1PSYDqlRmSVWF2vRpnU6YeaadBCOTgMX5xMOy/0lM47QdxZkOOw90hIobE0RTuX9WDIuDvWpTWeRdS7pUwsH40CK0nyesW/rFrs+zPKORc2HvCYFadzYqpKQxVUm7ECatRdPWyL055ACxrWSvaKmxsD1KJENt3x9n6fgmGqhNWoygi0GV6/n36oYMmIz0lnOiPELZdtrxUy/vj/dx63gnOLm/0U+MnKOgW3r2vQ7ZVSi5ttrn/gMEHR8uY0zbGPRf5Zxf8OwraYd43f414fVBw244KLmtyCZfYHOwsI5nbEuhUXvBuvMmErnej/Oj3DVx9QYoHl/cA0JkIbkYs+Dam7GfluD3AvW+W0NVZycVzE9SbkpwnWD87mDIp+DE3HQVTa2bDojAT84exPcn+IZdDw8F1kIJHHiZBLBcFzyWasUDVMRSHAx/prLxpq1xwk0V3Gh5c4F/WC67peORbYvtgmnX7kiQigks2l9NammXxtBwzrrNOXvfCh5t9q3FfbLEHttVxcE8td+7WaDt/mGROpT8tGMqCrsAPzs8ybitsOOSP/euL6HPaEpRXODRVTjpII1mNeETSaIIiJPHeAeIhh9604IEDhR93+FJne/sekzu7DbpTYI0b1MRyLKzzSNmCV3oFw91RBgYi7LrN48hoN7dsnsJfv3tp3rF9ES1T22TKAF2EeK4rwW/3lDP17DF+eAEMWWnW/66J8bSG6wkMbdLXs13Ba99bkld8X/qSHDiK6bggBJc0NdG6dBBVQG/aoyaiY+VUlAmJp4ItFCKGS1VlLu94vtS0MvUugWEg6s4i66a4+bEWDicMTM3DUCJsfeuXCCkxohJz1GbJoiGuX5zIO55vvXMos1nGvCg7BlOcU2HRPWEwO+ZSaTSQckZ5e9jijvYmcOD527vw8OhPK8yJX3/aGnwTvfPE49L2BIdGVOJhSZkB7yUFlgtLm3M4EtoqFDzpMmZPdhlNSFrLbjxtDb4dX+zu1TA1jytb02RdWNzg8f05NcyslGRdQaNpM2TlEELQVFLPdFlCi6XmFcu3TH/t0edlNqviuIK6hhRjr2qUfsll9bwJxmzB5c0ulc+/jKiPMzi3jRq9EaEvzSu+b5lWFElZWQ4z7DDQa5JqDbF0dgpdkcyLZ6jq64cV19Iz5yyq97+Tt2BfRSMhndKI1adpmjpObW2GbX+r4L2kSlO0EhFrYXfvCY6PhRm+cDZyNP+Z0df17oon/iDjLePUhmHnSzXEEhaZVgNd85g1Z4SfXGTSkxrild4I9abHspb8XsP3ddFbXp3BzSkMCY/Wcyfo6zGRQ5LIhMPQFIP1L+e4uDFCLEzegiGA0621O56S2bRGLqsyPBwmk1YImx5h3cNxBVVVWZbMHWN+TZa51affoyGAc0TbUgFBKOxSU5MhndJpmJ7EyaqMD4dBSHa+VcbAueN5xwhkD7fmma1S0SSq7qLqHmbII6JBwhI4lkrfh6W8sG7ZF6B7nELH8pViYtRA+dfoliNwPBCKRCiSr16c4NY9+XePQEQDbFjRi+sIVN1DVSUhBcp0mB5zKTMkc2P57wUC2zLXpRQem9IJUtI3rY5fvV1K5+EKGqYnqYlIbmjL378OTLQbb0GtbGbfSA+vHwthe+C5gmS/yV7H+s8DfAaBlYemXCoyWMwPhfiOdoyRUZ3EuIHrKpxf5xY0dmCiAUz9KkFJjPTG3Tg5la9/pZ9INIeuFLa3DdZIBoS+VNy9r0NaR3UGLYua6iz7BwrLVaCZ/pj75q8RmbRGcmIyR9HQGfTu4ZINL8pvdDx3BikuUuR/wz8BP5MUOG5tAsMAAAAASUVORK5CYII=",[[-33.941641403099,18.3003829808018],[-34.3371886093941,18.5267584324]],1,null,null,"lambda"]},{"method":"addRasterImage","args":["data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAC0AAABfCAYAAAB88E5QAAAFoUlEQVR4nO2a3WscZRTGf2f2I7ubz1ZiSAs1RbGtQrUifoEIpdJGaGsjLUiRgoJ/QL2st8WCF9IbRRRBcuFF0jZabdVGBBEUiuhNWyoFjWBr05C03WyyO7szx4uZ2Z2tTTEz7zS9mAc2mx2yc559ct7znvO8CylSLAlJ4qY6N6oU+7wAhZ3GY1imbzg+8KKGX9sf7NOl/jYqjJL+RLYqQGPiFJSv8Y97hfFtu0yGAAySnj6wQwGGhgpcGq3gnPme87MdfPbDKlMhmjBCWqc/UoBNW0p0loRSj4N7vcrm+1zccw7bDp02miKxSat9SgFKfQ3sukuuoEgGJJ/h7HSGck8Hk4eHjS7GbJw3q31aQaHrfqxSloH+LIObFik8swbZe4CXZGsi1SkWaQDEgmyOji0DrOMqxYOvIv2vC3xqgN4SIaO+0VMZEKH+4cdcnlxkaOKrRJS9FZFzWvLDAorkdojV00GtkjHJK0WKFClSpEiRIkWKFCkMQSvHjbtH/wfxLYSFiRUhHgtaHtOVUjyFzh9TrZxITH3jPoVWjiuqnomDgiqoi3TvNRYrGVO9qXKLNACZHKginXtixTVuqntQEPHUvg29uKmTCGnpHBHchq+wBVbGe6iZNI9vQC4FsQABy/LJegZr3NSAJElDKzWCRYkZpRN1ObU8rliWn9sWUtx1V1zVFClSpLgL0NqXqtWTiffWseqmVk4oIt5tJLT7la/5x3LJIHLvcdsxK7iSycVuiu6ESKSbhC3L7+YItaCC0zfAfK5hkGY7oikddGuqfksR9BXe9Yxdoztvrum/Fcsm3fZvD8hKiJ86xlrQpRCvyxNahCX0+SXZvmjZSkvnHmkjG6SHun5Oa+JKR8ppKe0RXMdP4WC0ksQVDhB93LoTwWCQTQjxNpeg9AX57KeHlOKPVHdCrMFWSi975MI5rYo6k4kmtYFpPORriMUVpr2r9a8TI25gsJW2JBvMrvMtg4QsFUwoHW6YRDzVnXqiZc+MhRAorUA2i2OBJYLrfqfiOnDzKrL6NWOLM7bSUtzdTsapk5EcYlcRp968fGl4pzHpzSSeaqsPcV1wGzT9DnWhu596TVg4NGyEuBHSbqGLVjOtUKu09yKqDD3vIJZQPbI7NnEjeXbTHlNLMhQzPVjVCqCQK3hVpGGD69DsSfyFG8dtMqJ0zirQcG1u2ldZzOGpnMnjCO1bemBKArr4eWTFjVSPgubpyJcQhXJjFrLduDhksP5b+sT/ESNJjJCW3HbRhQlFhG4pAnVqb7/H9SmLwaMj0DMQitjB/i8qzP2UjxzP4Lblq6cOqEt2bReuK8wdGkPP/QyNGqBQr7H90QVytsO+N8Yi6W3Qn9bQk5J96iHWjKxF+tZgvz9KrlqD7k7kwcd4uLfO3OoC6/+8ESmSsV2qzVJoDrrSqtWNGjp1gcaP58k+t5GZDesZvdjBW4/vXzYHc6SDgVdoLbLwNKOuV/UunOXKu78xeHAz1tNHIsU3SPq4tnV8SouwU/eO41zHq93laWTgzZWt0x4CiYOuL7iukPUrhZX1FC/2ovPjkYueMdLelB6SOXSihRsMCXhK5wrt2/wyYfZ0S0P9RyB1m5ET2sY7R+6F9PAPPcON0+2gCvPXcC8eXtltfEmo+inefgBKaRWSLUS+bUIH+se0rXqEv5Hgh4xzcpuYP9Esgd4r2qqLeL9H9UcSTI/wLhO+BucWZ2I5aAmSDlTFey72gl1Bfz/LhkeeJZfZcW9Uj3Z4Smuxl1mrTMOtwfwM1Y1PMFv7O9adEyPtbTYWUp1ndaafRafML47Lyakyl27EC5ug0r7X59TBXqDbtvh1Js+ZqQwNjbf+EyUNIF2vCKp8M1vBUXig1+XEH7lY90ycNHjEL1cs/poXZm2YvR591LrrePKdSX3h6LfptyVTpLgF/wLgRijPgQ+lqQAAAABJRU5ErkJggg==",[[-33.941641403099,18.3003829808018],[-34.3371886093941,18.5267584324]],1,null,null,"A"]},{"method":"addLayersControl","args":[["alpha","gamma","lambda","A"],[],{"collapsed":false,"autoZIndex":true,"position":"topright"}]}],"limits":{"lat":[-34.3371886093941,-33.941641403099],"lng":[18.3003829808018,18.5267584324]}},"evals":[],"jsHooks":[]}</script>